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ABSTRACT 

The present status of the solar neutrino problem is reviewed. The strongest 
motivation for new neutrino physics comes from the complete phenomenological 
failure of astrophysical solutions: (1) The standard solar model is excluded by 
Qh each of the solar neutrino experiments. (2) The combined results of Homestake 

\ and Kamiokande are incompatible with any astrophysical solution. (3) Even if 

If} ' the Homestake results were ignored entirely, there is no realistic solar model 

CN ■ so far that describes the Kamiokande and gallium data simultaneously within 

\ their experimental uncertainties. On the other hand, the MSW solutions provide 

' a complete description of the data, strongly hinting at the existence of mass 

s * ■ and mixing of neutrinos. The MSW predictions for the future experiments are 

, given in detail. It is stressed that solar model uncertainties are significant in 

■^j- | determining the oscillation parameters and in making predictions for the future 

experiments. Especially, a 8 B flux (or Sn) smaller by 20% than that of the 
q>^ . Bahcall-Pinsonneault model can make the interpretation of the charged to neutral 

\ current ratio measurement in SNO ambiguous. 

Oh" 

i 1. Introduction 

There is a mounting evidence that the solar neutrino flux is considerably less than 
the theoretical predictions of the standard solar model (SSM). The SSM predictions by 
Bahcall and PinsonneaultS and by Turck-Chieze and LopesS are listed in Table [I] along 
with the experimental results. The recent SSM by Bahcall and Pinsonneault uses the 
updated input parameters and includes the helium diffusion effect. The observed rates 
of the Homestake chlorine experiment sBi Kamiokande water Cerenkov experiment ,@i 
and the combined result of the SAGeI! and GALLExi gallium experiments are 3.1, 
5.5, and 4.6 standard deviations smaller than their predictions, including the theoretical 
uncertainties. The deficit is 1.3, 2.9, and 3.7 <r, respectively, when the data are compared 
to the Turck-Chieze-Lopes SSM. The differences of the two calculations are understood 
as differences of input parameters, and also the Turck-Chieze-Lopes model does not 
include the helium diffusion effect and assigns more conservative uncertainties. 

The solar neutrino data are graphically compared to the SSM predictions in Fig. [I]. 
The solar neutrino problem is summarized in two aspects: 

*Based on a talk given at Solar Modeling Workshop, Institute of Nuclear Theory, University of Wash- 
ington, Seattle, March 24, 1994. 
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Table 1. The standard solar model nredictions of Bahcall and Pinsonneault (BP SSM)EI and of 
Turck-Chieze and Lopes (TL SSM),@ along with the results of the solar neutrino experiments. 



BP SSM TL SSM Experiments 



Kamiokande" 


1 ± 0.14 


0.77 ± 0.19 


(0.51 ± 0.07) x 5.69E6 /cm 2 /sec (BP SSM) 


Homestake 6 (CI) 


8 ± 1 SNU 


6.4 ± 1.4 SNU 


2.32 ± 0.26 SNU (0.29 ± 0.03 BP SSM) 


SAGE C & GALLEX d (Ga) 


131.5 tl SNU 


122.5 ± 7 SNU 


78 ± 10 SNU (0.59 ± 0.08 BP SSM) 



° The combined result of Kamiokande II [ 0.47 ± 0.05 (stat) ± 0.06J>ys) BP SSM, 1040 days]! 
and Kamiokande III [ 0.55 ± 0.06 (stat) ± 0.06 (sys) BP SSM, 627 days]0 is 0.51 ± 0.04 (stat) ± 0.06 
(sys) BP SSME The 8 B flux in the Bahcall-Pinsonneault (BP) SSM with the helium diffusion effect 
is 5.69xl0 6 /cm 2 /sec. 

6 The result through June, 1992 (Run 18 - 124) is 2.32 ± 0.16 (stat) ± 0.21 (sys) SNU B 

c The preliminary result of SAGE I (through May, 1992) is 74 ± 19 (stat) ± 10 (sys) SNU.H 

d The combined result of GALLEX I and II (30 runs, through October, 1993) is 79 ± 10 (stat) ± 6 

(sys) SNU .B 



• Assuming the standard properties of neutrinos, the SSM is excluded. Neither the 
Bahcall-Pinsonneault or the Turck-Chieze-Lopes model (or any other SSM) are 
consistent with the observations. 

• More importantly, the smaller observed rate in the Homestake data relative to 
Kamiokande is inconsistent with astrophysical solutions in general. When the 
neutrino flux is reduced by changing the solar model, the 8 B flux is generally 
reduced more than the 7 Be flux, and therefore the Kamiokande rate should be 
suppressed more than the Homestake rate, contrary to the data. 

From a particle physicist's point of view the most exciting interpretation of the 
problem is that the solar neutrino deficit is a manifestation of nonstandard neutrino 
properties, such as mass and mixing assumed in the Mikheyev-Smirnov-Wolfenstein 
(MSW) mechanism. This possibility is, however, one of three generic possibilities, and 
any combination of the three is possible. 

• Experimental Solutions: One or more of the experiments is wrong. 

• Astrophysical Solutions: The SSM is wrong. Nonstandard input parameters and/or 
nonstandard processes that are not included in the SSM will explain the deficit 
of the solar neutrino flux. 

• Particle Physics Solutions: The standard model of the electroweak theory is 
wrong. Nonstandard neutrino properties are causing the deficit. Proposed so- 
lutions include the MSW effect, vacuum oscillations, a large neutrino magnetic 
moment, neutrino decays, violation of the equivalence principle, etc. 
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Fig. 1. A graphical summary of the solar neutrino problem. The experimental rates of the 
combined SAGE and GALLEX, Homestake, and Kamiokande experiments are shown with 
respect to (a) the Bahcall-Pinsonneault SSM and (b) the Turck-Chieze-Lopes SSM. The data 
points are plotted according to the energy sensitivities of the experiments, but only crudely. 
The theoretical uncertainties are also displayed as error bars (1<t). 

In this paper I will focus my discussion on the theoretical aspects of the solar neutrino 
problem. For astrophysicists what is at stake is the understanding of the Sun, which is 
the best measured main sequence star. For particle physicists, it is new neutrino physics 
such as neutrino mass and mixing. These neutrino properties require a (long-waited) 
extension of the standard model and provide more information on the fermion mass 
spectrum and fermion mixing, whose theoretical understanding is totally lacking. 

It is not my intention here to discuss the reliability of the experiments explicitly 
other than mentioning that no defect has been found so far in any of the existing 
experiments. Instead I will use the experimental data as given and will quantitatively 
test various theoretical ideas. I also address my objection to the theoretical analyses 
in the literature done by selecting a subset of available data set (e.g., discarding the 
earlier Homestake data). This kind of data selection is unfounded and misleading, and 
often defocuses the issue we are facing. 

This paper is organized as the following. In Section |2], astrophysical solutions are 
discussed for the cool Sun model and for the model independent analysis. In fact the 
most compelling reason for considering particle physics solutions comes from the com- 
plete failure of astrophysical solutions to describe the data. It will be shown that the 
observed lower Homestake rate relative to the Kamiokande rate essentially excludes 
all astrophysical solutions. The lack of a viable astrophysical explanation motivates us 
to consider particle physics solutions. In Section H, the MSW solutions are discussed 
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as the simplest particle physics solution. The results of the updated MSW analysis 
are presented. The analysis includes the Earth effect, the Kamiokande II day-night 
data, the SSM theoretical uncertainties and their correlations, a proper joint analysis 
of the experiments, and an improved determination of the confidence levels. Assum- 
ing the SSM, the experiments allow two parameter regions at 95% C.L.: one in the 
small-angle (nonadiabatic) and one in the large-angle region, the former solution giv- 
ing a better fit. There is also a small-angle (nonadiabatic) solution for oscillations to 
sterile neutrinos, but none in the large-angle region even at 99% CL. The MSW effect 
can also be considered in nonstandard solar models. The allowed parameter space is 
obtained for nonstandard core temperatures, nonstandard 8 B fluxes, and various ex- 
plicitly constructed nonstandard solar models. It is stressed that, once the MSW effect 
is established in the next generation solar neutrino experiments, the uncertainties in 
the neutrino parameters will be dominated by the solar model uncertainties, and the 
calibration of the neutrino flux in the presence of the MSW effect will be a major 
theoretical issue. Assuming the SSM and the existing data are correct, the MSW effect 
provides robust predictions for the next-generation experiments such as SNO, Super- 
Kamiokande, BOREXINO, and ICARUS, and those predictions are given in detail in 
Section [| I argue that the prediction for the charged to neutral current ratio in SNO 
is solar model dependent, and a correct estimation of the 8 B flux uncertainty is crucial. 

2. Astrophysical Solutions 

The SSM calculations have been refined over the years, and the recent model 
by Bahcall and Pinsonneault improved on earlier calculations by incorporating the 
new Livermore opacity and updated element abundances, and also by including the 
helium diffusion effect for the first time in the SSM calculation! The model is in good 
agreement with the helioseismology observations. They also showed that the various 
SSM calculations give essentially the same flux predictions if the same input parameters 
were used, proving the reliability of the calculations within the SSM framework. More 
recent attempts to incorporate mixing processes, heavy element diffusion, and rotations 
were reported by other speakers at this workshop. Erllil 

Despite those successes and improvements of the SSM, one has to depart from 
the SSM in order to search for an astrophysical solution of the solar neutrino problem. 
It is well known that the neutrino flux prediction responds sensitively to the changes 
in some of the solar model input parameters. For example, opacity or heavy element 
abundance smaller than the standard values can achieve lower neutrino fluxes. Astro- 
physical processes that are not included in the SSM calculations, such as rotations, 
core magnetic fields, or gravitational settling of elements, might significantly reduce 
the neutrino flux. Those effects usually work by reducing the core temperature, and 
the neutrino fluxes are sensitive to the change in the core temperature. 

Low energy nuclear reaction cross sections are another source of flux uncertainties. 
The low energy cross sections important in the SSM such as Sir, S33, S34, and Sn have 
never been directly measured at the low energy equivalent to the solar temperature, and 
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there could still be unquantified theoretical uncertainties in extrapolating the measured 
cross sections to the lower energies. Such uncertainties might be correlated among 
different cross sections. 

In particular Sn, to which the 8 B flux magnitude is directly proportional, is a 
poorly measured quantity and I believe that its quoted uncertainty is underestimated. 
The standard value (617(0) = 22.4 ±2.1 eV b) recommended in Ref. Hl5l , ll6]l is the 
combined result of six measurements, but is controlled mostly by the two inconsistent 



measurements with relatively small uncertainties: Sn(0) = 25.2 ± 2.4 eV b of Ref. |17 
and 20.2 ±2.3 eV b of Ref. |18j. From the differential cross section plot of the two data 
shown in Ref. ||15|| , it is clear that the two measurements are inconsistent: if one result 



is correct, the other must be wrong. The true value can be 20.2 — 2.3 = 17.9 eV b or 
25.2 ±2.4 = 27.6 eV b. My objection for the combined value is that its uncertainty does 
not cover the entire possible range. As I will show in Section [2.2|, the lower choice of Sn 



does not solve the solar neutrino problem. However, underestimating the uncertainty 
does change the obtained MSW parameters, and therefore change the MSW predictions 
for the future experiments. The charged to neutral current ratio experiment in SNO is 
the gold-plated measurement to establish neutrino oscillations, but, given the currently 
projected systematic uncertainty of 20%, a 8 B (or Sn) lower by 20% than the Bahcall- 
Pinsonneault value can make the conclusion of the measurement ambiguous. 

The worst scenario is that those astrophysical effects and nuclear physics effects 
might be taking place at the same time. In fact changes in nuclear reaction cross 
sections such as Sn can cause a change in the core temperature. This complexity and 
uncertainty of the problem allow theorists to come up with new astrophysical solutions 
now and then, such as a recent claim of the standard solar model solution to the solar 
neutrino problem.lll (Fatal flaws of the claim are discussed elsewhere .&H'@) 

Despite all those worrisome aspects of the solar model calculation, it can be shown 
that the data cannot be explained by astrophysical and nuclear effects. In the following, 
I will discuss two analyses of nonstandard solar models: the cool sun model, follow- 
ing the analysis of Ref. [^3|,0] (see also Ref. p5|,|9|, |26|j27|] ), and the model independent 



analysis, following Ref. |p8fl , which depends essentially only on the data and provides 
powerful conclusions. The cooler Sun analysis shows that the data exclude a wide class 
of nonstandard solar models characterized by lower core temperatures; the model in- 
dependent analysis shows that, under minimal assumptions, the data are incompatible 
with any astrophysical/nuclear solutions. 

2.1. Cool Sun Model 

A class of nonstandard solar models can be characterized by lower core tem- 
peratures.il@ The validity of this general description was explicitly demonstrated by 
Degl'Innocenti for solar models with large Sn, low opacity, low Z, and low age.0 The 
neutrino fluxes are described by thepower laws of the central temperature (Tc) ob- 
tained from 1000 Monte-Carlo SSMsBfl: 

<f>(pp) ~ Tc 1 - 2 , 0(Be) ~ T c , and 0(B) ~ T c 8 , (1) 
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Note that these relations are obtained within the SSM calculation, and there are uncer- 
tainties in the exponents in the power law, especially when Tq is changed beyond the 
SSM uncertainty (Tc — 1 ±0.006). (The exponents obtained by Degl'Innocenti is —0.6, 
9, and 21 for the pp, 7 Be, and 8 B flux, respectively.0) Especially, the power law for the 
pp flux is only valid within the SSM uncertainty, and a more reliable Tc dependence is 
obtained by the luminosity constraint, which determines the total number of neutrinos 
generated by nuclear reactions. Taking account the energies carried off by neutrinos, 
the luminosity constraint yields 

<f>(pp) + 0.979 0(Be) + 0.955 0(CNO) = 6.57 x 10 10 cni'V 1 , (2) 

where the fluxes are described as functions of Tc- We use the power law (Eqn. |l|) for the 
7 Be and 8 B flux and take the exponent of the CNO flux as 22.0 4>{pp) is then obtained 
from Eqn. |2|. The central temperature is in units of the Bahcall Pinsonneault SSM value 
(Tc = 1 = 15.57 x 10 6 K). Using the Tc dependence, one can write the expected rates 
of the Kamiokande (i^Kam), chlorine (-Rci), and gallium (i?Ga) experiments relative to 
the SSM predictions as functions of Tc'- 

R Kam = (1 ± 0.089) T c 8 (3) 
ifci = (1± 0.033) [0.775 (1 ± 0.089) 7£ 8 + 0.150(1 ±0.036) T 8 + small terms] 

R Ga = (1 ± 0.04) [ 0.538 , + 0.272 (1 ± 0.036) T c (4) 

tPWP) SSM 

+0.105 (1 ± 0.089) T c 8 + small terms ]. 

Our analysis includes the uncertainties (shown in parentheses) and their correlations 
from the nuclear reaction cross sections and the detector cross sections. 

When the data are fit to Tc separately, the Kamiokande and Homestake data 
require Tc = 0.96 ± 0.01 and 0.92 ± 0.01, respectively. The gallium rate of 78 ± 10 
SNU requires T c = 0.71 ± 0.14, whose central value is unrealistically low compared to 
the SSM value (T c = 1± 0.006). Moreover, the temperatures obtained are inconsistent 
between different experiments. When all data are fit simultaneously, the best value is 
Tc = 0.93 ± 0.01, but only with an extremely poor fit; x 2 — 15.7 for 2 degrees of 
freedom, and statistically this model is excluded at 99.96% C.L. The results of Tc fit 
to various combinations of the data are listed in Table |2], along with the expected rate 
for each experiment from the obtained temperatures. It should also be noted that there 
is no consistent temperature even if we ignore one of the three experiments. 

We have reanalyzed the data with various other combinations of the exponents of 
the Tc power laws, and have found that the as long as the 8 B flux is more temperature 
sensitive than the 7 Be flux, the combined observations cannot be described by lower 
core temperatures. The fit has been repeated with all theoretical uncertainties tripled. 
The result is somewhat modified, but the conclusion remains the same: the cool suns 
do not explain the solar neutrino observations. 
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Table 2. The results of the Tc fits. The expected rates for the three experiments are listed 
for each obtained Tc] the Kamiokande rate is in units of the Bahcall-Pinsonneault SSM and 
the chlorine and gallium rates are in units of SNU. 





T c 


± 


AT C 


x 2 


C.L. (%) 


Kam 


CI 


Ga 


Kam 


0.96 


± 


0.01 







0.51 


4.2 


111 


CI 


0.92 


± 


0.01 







0.22 


2.3 


99 


Ga 


0.71 


± 


0.14 


0.08 




0.002 


0.16 


81 


Kam+Cl 


0.93 


± 


0.01 


10.6 


99.98 


0.29 


2.7 


102 


Kam+Ga 


0.96 


± 


0.01 


9.0 


99.7 


0.44 


3.8 


109 


Cl+Ga 


0.92 


± 


0.01 


4.2 


96 


0.23 


2.2 


99 


Kam+Cl+Ga 


0.93 


± 


0.01 


15.7 


99.96 


0.28 


2.6 


102 



2.2. Model Independent Analysis 

A more powerful conclusion can be obtained just from the observations, without 



referring to any specific nonstandard solar models. Here I follow the analysis of Ref. p8 



but a similar approach can be found in Ref. [32 . 

Suppose we do not know anything about the theory of the Sun; what will the 
solar neutrino data tell us about the star? We have formulated the question by using 
four relevant neutrino fluxes (pp, 7 Be, 8 B, and CNO) as free parameters fit to the ob- 
servations. Except for the luminosity constraint, there are no relations imposed among 
the fluxes as in the cool sun model. It is not my intention to advocate such a model as a 
realistic description of the Sun or to claim that it is consistent with other observables, 
but to illuminate the problems of the astrophysical solutions when they are confronted 
by the observations. 

While using the four fluxes as free parameters, we made three explicit assump- 
tions: 

1. The Sun is in a quasi-static state, and the solar luminosity is generated by the 
ordinary nuclear fusions. That is, the Sun shines as it constantly consumes its 
nuclear fuel through the pp and CNO chains. This imposes a relation among the 
fluxes (Eqn. |). 

2. Astrophysical mechanisms do not distort the neutrino energy spectrum at the 
observable level. We are allowed to change, for example, the amplitude of the 
8 B flux freely, but not to deplete only the lower energy part of the spectrum 
in order to reconcile the Homestake and Kamiokande rates. It was shown by 
Bahcall that possible distortions of the spectrum due to such astrophysical effects 
as gravitational red-shifts and thermal fluctuations are completely negligible.0 
Note that particle physics solutions are often energy dependent and do distort 
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the spectrum: the MSW small-angle (nonadiabatic) solution is a typical example. 



3. The experiments are correct, and so are the detector cross section calculations. 



The results of the analysis are presented in the 2-dimensional plane of the 7 Be 
and 8 B flux in units of the Bahcall-Pinsonneault SSM values. The analysis includes 
the experimental uncertainties and the uncertainties from the detector cross sections 
and from the minor fluxes. When the data are fit, the pp and CNO fluxes are changed 
freely for each 0(Be) and 0(B), subject only to the luminosity constraint. This plane 
shown in Figures |2] and |3] essentially represents every possible astrophysical solution: 
from the minimum rate model (0(Be) = 0(B) = 0) to the SSM (0(Be) = 0(B) = 1). 
That is, any claims — in the past, present, and future — of the astrophysical solutions 
are represented in Figures and |3] as long as they satisfy the minimal assumptions. 
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Fig. 2. The allowed fluxes when the Homestake, Kamiokande, and the gallium results are fit 
separately. This is an updated analysis of Ref. [28|. 



When each experiment is fit separately, the Kamiokande data constrain the 8 B 
flux, and the Homestake and gallium experiments constrain both the 7 Be and 8 B fluxes; 
constraints at the 90% C.L. from each experiment are shown in Fig. ^| In this model 
0(Be) = 0(B) = is completely consistent with the Homestake and gallium observa- 
tion, since we can always adjust the pp and CNO fluxes to the observed rates. As seen 
in Fig. |2|, each of the experiments is inconsistent with the SSM prediction. 
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Fig. 3. The fluxes allowed by the combined data. Also shown are various standard and 



nonstandard solar models. This is an updated analysis of Ref. [28 |. 



When the observations are fit combined, astrophysical solutions in general are 
excluded, and the argument comes in three steps. 

The x 2 value of the joint fit is very poor. In fact if (nonphysical) negative fluxes 



were allowed, the x minimum is obtained for 0(Be)/0(Be" 



SSM 



-0.5. Within the 



physical region, the best fit is obtained for 0(Be)/0(Be)ssM < 0.08 and 0(B)/0(B)ssm = 
0.37 ± 0.04 (la), but only with x 2 — 4.1 for 1 d.f. Statistically this fit is excluded at 
96% C.L. In other words, any astrophysical solutions consistent with our minimal 
assumptions are excluded at 96% C.L. by the observations. 

Secondly, even if I allow the 4% chance and accept the solution, the obtained 
fluxes do not physically make sense at all. Since 8 B is produced by the reaction p + 
7 Be — > 8 B + 7, essentially any reduction in the 7 Be flux should cause at least an equal 
reduction in the 8 B flux, and additional temperature sensitivity reduces the 8 B flux 
more. Therefore, unless the uncertainty of the 7 Be electron capture rate is grossly 
underestimated or there is some independent mechanism to suppress only the 7 Be 
neutrinos, any realistic model should be in the region 0(Be)/0(Be)ssM > 0(B)/0(B) SSM 
in Fig. H contrary to the observations. 

Finally various standard and nonstandard solar models are plotted in Fig. |3|. The 
Bahcall-Pinsonneault SSM with 90% C.L. uncertainty^ the Bahcall-Ulrich 1000 Monte 
Carlo SSMs@ the Turck-Chieze-Lopes SSM,i the solar models with opacity reduced 
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by 10 and 20%£ 4 ] the large Su models by Castellani et al.fB the low Z model by Bahcall 
and Ulrich,E3'E3 the WIMP model by Gilliland et al., 00 and the Dar-Shaviv model.0 
Also shown are the_power law dependence of the fluxes to T c and Su obtained by 
Bahcall and Ulrich.0'0 None of those models is even close to the observations. 

This complete failure of our model forces us to conclude that at least one of the 
original assumptions is wrong. We have repeated the analysis without the luminosity 
constraint; (j the exclusion of the model is somewhat weaken (the exclusion is at 95% 
C.L.), but the conclusion remains essentially the same. Therefore we have to conclude 
that either of the other two assumptions is incorrect: (a) there is a distortion of the 8 B 
flux, that is, new neutrino physics, or (b) one or more of the experiments is wrong: in 
particular, either the Homestake or Kamiokande data must be grossly wrong. [[ 

2.3. What if Homestake data were ignored? 

As discussed in the previous section, the most compelling motivation for the 
particle physics solutions such as the MSW effect comes from the comparison of the 
Homestake and Kamiokande data, and, if one of the two were ignored, the astrophysi- 
cal solutions cannot be dismissed outright. Some argue that we would lose motivations 
for new neutrino physics if the Homestake data were ignored. This, however, is mis- 
leading: at present we have no realistic solar model that are consistent both with the 
Kamiokande and gallium results within their uncertainties. 

I address the issue by considering the solar models with the 8 B flux consistent 
with the observed Kamiokande rate.0 The Monte Carlo SSMs with such a constraint 



were discussed in Ref. ||25|| , but I also allow nonstandard solar models. The neutrino flux 
varies significantly from model to model in nonstandard solar models in general, but the 
fact that about half of the SSM 8 B flux is seen strongly constrains the nonstandard solar 
models, and yields predictions for the gallium rate and the helioseismology observations. 



* The only motivation for this is to allow for the possibility of models with extremely rapid variation 
in the solar core (i.e., faster than the 10 4 year time scale for photons to diffuse to the surface), such 
as the model of Ref. |J7| . 

t It was argued by Davis@0 that the 8 B flux observed in Homestake and Kamiokande is consistent, 
and it is explained by the mixing model, which assumes an ad hoc mixture in the solar interior.H 
The comparison is, however, made with highly selective data sets: the data only during 1986 - 1990 
are used. Furthermore the Kamiokande II rate is taken from the data set with the energy threshold 
9.3 MeV, ignoring the data with the threshold 7.5 MeV. The data used in the comparison is about 
15% of the total Homestake data and 30% of the currently available data from Kamiokande II and 
III. Given the consistency among the Kamiokande results with different thresholds, Em this selection 
of the Kamiokande data is unjustified. The selection of the data during a particular time period is 
justified only when a time variation in the flux is assumed. The mixing model, to which the selected 
data were compared, does not predict time variations, and therefore the theory should be compared 
to the averaged data. The mixing model (with a slow mixing for 40% mass) predicts the Kamiokande, 
Homestake, and gallium rates to be 0.44 BP-SSM, 3.6 SNU, and 104 SNU , while the averaged data 
(without any selection) are 0.51 ± 0.07 BP-SSM, 2.32 ± 0.26 SNU, and 78 ± 10 SNILjespectively. 
This mixing model is inconsistent with the observations. As discussed by Merryfield,cil the mixing 
model is also disfavored by the helioseismology data. 
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I list in Fig. f| the expected gallium rates from various nonstandard models that are 
consistent with or close to the 8 B flux observed by Kamiokande: the model in which 
S17 is adjusted from the Bahcall-Pinsonneault model to the observed Kamiokande rate 
(0.51 ± 0.07 SSM), the low Tc model, the low Tq model with Sn increased by 30%, 
the model with S34 reduced to the Kamiokande rate, the Dar-Shaviv model,El the low 
opacity model (— 20%), the model with Sn increased by 30%il and 25%,! and the 
mixing models by Sienkiewicz et al.0 The uncertainties include the uncertainty of the 
8 B flux due to the Kamiokande uncertainty (14%), but the dominant contribution is 
from the gallium cross section uncertainty. From this list, I obtain 

Gallium rate consistent with Kamiokande > 100 SNU, (5) 

while the combined gallium rate of SAGE and GALLEX is 78 ± 10 SNU. 

Obviously it is still premature to rule out those nonstandard solar models given 
the current large experimental uncertainties. But the observations and the theoreti- 
cal predictions are showing discrepancies. Also other solar observations might exclude 
those models. The helioseismology observations provide information on the sound speed 
throughout the radiative and convective zone (but not in the core region where the the- 
oretical and experimental uncertainties are still large). As discussed by other speakers 
at this conference, the observations give significant constraints in solar model building: 
the SSM is consistent with the helioseismology data at the 1% level, and any large 
departures from the SSM will be likely to destroy this agreement. In fact, some of the 
nonstandard solar models listed in Fig. [1 and also other nonstandard models such as 
the low Tq model,@ the large Sn model]! the mixing model,0 and the low Y modefi 
are in conflict with the helioseismology observations. 

While waiting for the gallium detectors to be calibrated and the statistics of the 
gallium data to improve, this line of theoretical investigations should be encouraged. 
I hope that the possibility of solar models consistent with the Kamiokande and the 
gallium data will be pursued, and the list of nonstandard models as in Fig. |3| will be 
extended to obtain the minimum rate expected for the gallium experiments. Also, each 
and every model should be confronted with the helioseismology observations. If one 
could successfully find a model consistent with the Kamiokande and gallium as well as 
the helioseismology observations, it could be the solution to the solar neutrino problem: 
the SSM and the Homestake experiment might be wrong after all. If, however, one could 
not find such a model after entirely discarding the Homestake data, our motivation for 
considering particle physics solutions will be complete. 

3. MSW Solutions 

There are many proposed particle physics solutions to the solar neutrino prob- 
lem: the MSW effect with two or more flavors, vacuum oscillations, a large neutrino 
magnetic moment, neutrino decays, violations of the equivalence principle, etc. I dis- 
cuss the two state MSW oscillations,!^! since it requires only the minimal extension of 
the standard model (mass and mixing of neutrinos) without contradicting with any 
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Fig. 4. The rate of the gallium experiments, the gallium rates of the SSMs, and of the 
nonstandard solar models which predict the 8 B flux consistent with or close to the flux 
observed in Kamiokande. See text for details. The nonstandard solar models consistent with 
Kamiokande predict the gallium rate i?Ga > 100 SNU. The MSW solution obtained from the 
combined Kamiokande and Homestake data predicts Rc a < 100 SNU, consistent with the 
data. 
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other existing observational constraints, and also without requiring any unnatural fine- 
tunings of the parameters. Furthermore, many extensions of the standard electroweak 
model with the see-saw mechanismB suggest neutrino mass consistent with the MSW 
parameter range,0 providing a theoretical motivation to look for neutrino mass in the 
solar neutrino experiments. I also consider the MSW oscillations for sterile neutrinos 
and the MSW solutions in the context of nonstandard solar models. The MSW analysis 
of the solar neutrino data was carried out by many authors,&H but here I follow the 



recent analysis of Ref . [^j,|>|,[Jl[] , which includes the theoretical uncertainties and their 
correlations, the Earth effect, the Kamiokande II day- night data, and the improved 
definition of the confidence levels. 

3.1. Earth Effect 

One interesting, but not-well publicized, aspect of the MSW mechanism is the 
Earth effedS on the solar neutrinos. The electron neutrinos produced in the Sun trans- 
form to another species of neutrino by the MSW effect in the Sun, but, for a certain 
parameter range, they can transform back to ^ e 's during the night due to the mat- 
ter effect as the neutrinos propagate through the Earth. Such an effect can be seen 
as a difference in the observed rate between day and night, or, more drastically, as 
signal variations as a function of time during the night since the oscillation proba- 
bilities depend sensitively on the neutrino path length in the Earth and also on the 
densities that the neutrino go through (i.e., the mantle only or both the mantle and 
the core). The Earth effect is only relevant for a limited region of Am 2 and sin 2 29, 
but, if measured, it would be an indisputable evidence for the MSW solution. Given 
the high statistics in the next generation experiments, the MSW parameters (Am 2 and 
sin 2 29) could be determined precisely by measuring time and energy dependence. Since 
the measurement depends only on the time (and energy) dependence, but not on the 
absolute neutrino flux, the determination of Am 2 and sin 2 29 would be free from the 
solar model uncertainties and also from the time-independent systematic uncertainties 
in the experiments; the uncertainties from the Earth density profile is negligibly small. 

The MSW transition is most significant when the resonance condition is satisfied: 

/\ ???, ^ 

— — cos 2^ = V2G F n, (6) 
2E 

where E is the neutrino energy, Gf is the Fermi constant; n = n e for oscillations to 
or v T , and n = n e — n n /2 for oscillations to sterile neutrinos, where n e and n n are 
the electron and neutron densities. The density profile of the Earth is approximately 
described as two step functions: the core (10 - 12.5 g/cm 3 ) and the mantle (4 - 5 g/ cm 3 ). 
A more detailed density profile can be found in Ref. ||56|| . For solar neutrinos energies, 
the parameter space relevant for the Earth effect is Am 2 ~ 1CT 7 — 1CT 5 eV 2 and 
sin 2 29 > 0.01. To compare the MSW predictions to the averaged rates of observations, 
one needs to average the v e survival probability over the night, since it depends on 
the path length and also on the density profile through which the neutrinos propagate. 
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The survival probability also vary each night due to the obliquity of the Earth rotation 
axis and should be averaged over a year. Also it differs for each experiment due to the 
detector latitudes. This averaging process requires intensive computations. 

The contours of the signal-to-SSM ratio for the Kamiokande experiment are shown 
in Fig. [5| for the day-time rate (a) and the night-time rate (b). For a certain parameter 
space, the night-time enhancement is huge: for example, the night-time rate is almost 
three times the day-time rate for Am 2 ~ 3 x 10 -6 eV 2 and sin 2 26 ~ 0.2. 

Real-time counting experiments such as Kamiokande are capable of observing 
the day-night effect. In fact the Kamiokande II collaboration has published the day- 
night data with 6 bins (one day-time data point and 5 night-time data points for 
different angles between the Sun direction and the nadir at the detector; the binning 
corresponds to selecting neutrinos with different path lengths in the Earth). The data 
show no enhancement of the signal at night, albeit still large statistical uncertainties, 
and exclude the parameter space shown in Fig. [5] (c). This exclusion of parameter space 
is insensitive to the solar model uncertainties and independent of solar models. 

The Earth effect has significant phenomenological implications to the existing and 
future direct-counting experiments. The day-night data are available so far only from 
the Kamiokande II experiment. The additional data from Kamiokande III is expected 
to double the statistics, and the updated analysis should further constrain the param- 
eter space. In particular, a significant day-night effect is expected in the large-angle 
solution of the combined observations, which might be observable if the Kamiokande 
II data and the (still unpublished) Kamiokande III day-night data are combined. With 
large counting rates in the future solar neutrino experiments, one should be able to 
observe a day-night effect at a few per cent level, and the parameter space in which 
the contours in Fig. [5] (b) are distorted from Fig. |] (a) should produce observable 
signals. If the day-night effect is observed, its time and energy dependence provides 
a clean determination of Am 2 and sin 2 26 without referring to solar model calcula- 
tions. The detailed predictions of the day-night effect for SNO, Super-Kamiokande, 
and Kamiokande obtained from the MSW solutions of the combined observations will 
be given in Section |3~6| . 

3.2. MSW Analysis 

Given a solar model, either standard or nonstandard, the MSW mechanism is 
a solid, calculable theory. It can be quantitatively tested and constrained, and the 
following is the questions to be addressed in the MSW analysis: 

1. Is the MSW hypothesis acceptable by the observations? 

2. If acceptable, what is the parameter space allowed by the observations? 

3. What are the predictions for the future experiments from the obtained parame- 
ters? 
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Fig. 5. The contours of the signal-to-SSM ratio in Kamiokande for (a) day time and (b) 
night time, (c) Also shown is the region excluded by the Kamiokande II day-night data with 
6 time bins. 
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The first question tests the whole idea of the MSW theory, including the assumption of 
neutrino mass and mixing. The answer to the second question will provide important 
information for the neutrino mass spectrum and mixing matrix of the lepton sector. 
And most importantly, the MSW theory is verifiable or falsifiable by the next generation 
experiments. 

To answer those questions quantitatively, it is important to include relevant the- 
oretical uncertainties in the analysis. Especially the 8 B flux uncertainty in the Bahcall- 
Pinsonneault SSM is 14% and comparable to the experimental uncertainties of Home- 
stake (11%) and Kamiokande (14%). In the Turck-Chieze -Lopes SSM, the 8 B flux 
uncertainty is 25% and dominates the experimental uncertainties. 

Equally important, but often ignored, are the correlations among the theoretical 
uncertainties. The leading uncertainty from the 8 B flux is strongly correlated among the 
experiments, especially between Homestake and Kamiokande. Also the uncertainties are 
correlated between different flux components. If, for example, the opacity were lower 
(or the core temperature were lower), both the Be flux and 8 B flux would be reduced. 
Therefore a careful estimation of the theoretical uncertainties and their correlations are 
essential in determining the neutrino parameters. In fact if the theoretical uncertainties 
were ignored, there would be no allowed solution in the large-angle region at 90% C.L. 
If the correlations were ignored, a significantly larger region would be obtained for the 
large-angle solutionis 

We have parameterized the SSM flux uncertainties by the uncertainties in the cen- 
tral temperature and in the relevant nuclear reaction cross sections (Su and S^). 113130 
It was shown explicitly that our treatment reproduces the flux uncertainties and corre- 
lations obtained from the Bahcall-Ulrich 1000 Monte Carlo SSMs.0 The parameteriza- 
tion method can be easily generalized to other solar models, for which the Monte Carlo 
estimations are unavailable. The 7 Be and 8 B flux distributions of the Bahcall-Ulrich 
1000 Monte Carlo SSMs and our parametrized uncertainty estimation are displayed 
in Fig. |]. The amplitudes and the correlation of the flux uncertainties are consistent 
between the two methods. 

One ambiguity in determining the neutrino parameters is the definition of confi- 
dence levels. The general prescription in this case is that, once the MSW mechanism 
is accepted as a viable solution, the probability density of the parameter space allowed 
by the observations is distributed throughout the log Am 2 — log sin 2 29 plane, and is 
given as P = iVexp(— x 2 /2), where N is the normalization factor and the \ 2 is calcu- 
lated for each combination of Am 2 and sin 2 2^.0 When the probability distribution is 
approximated as a gaussian around the x 2 minimum, the allowed region is determined 
by x 2 (Am 2 ,sin 2 2fl) < xLn + A X 2 with A X 2 = 4.6, 6.0, and 9.2 for 90, 95, and 99% 
C.L., respectively. The A% 2 is determined by the fact that the probability density is 
distributed in the 2-dimensional parameter space, and does not depend on the number 
of data points which enter the \ 2 calculation. The assumption of one gaussian distri- 
bution is questionable for the MSW analysis, which yields multiple x 2 minima in the 
parameter space. We have generalized the A% 2 method in the presence of multiple \ 2 
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Fig. 6. The 7 Be and 8 B flux distributions of the Bahcall-Ulrich 1000 Monte Carlo SSMs 
(dots). Also shown (solid line) are the flux uncertainties (at 90% C.L.) parametrized by 
the central temperature and the relevant nuclear cross sections. The two methods provide 
consistent flux uncertainty estimations. Taken from Ref. 



minima,^ and will use the improved definition in the joint analysis. 

We use the Petcov formulaElS to calculate the v e survival probability, and have 
compared the results to other analytic approximations (the Parke formula with the 
Landau-Zener approximation!! and the Pizzochero formula@), and found that the dif- 
ference in the combined fit is completely negligible. Also our calculations were compared 
to those of Fiorentini et al.,H who solve the neutrino propagation equation numerically: 
the agreement was excellent. 

Our MSW calculations include the MSW double-crossing, the neutrino produc- 
tion distributions in the Sun,0 the electron density distribution^ the detector cross 
sections,@iy and the energy resolution and the detector efficiency for Kamiokande J 

The uncertainties from the neutrino production distributions and from the elec- 
tron density distributions are studied by comparing three different solar models (but us- 
ing the same flux magnitude)0: the Bahcall-Ulrich SSM and the Bahcall-Pinsonneault 
SSM with and without the helium diffusion effect. No difference was observed. We have 
also changed by ±10% each the peak location of the production distribution, the elec- 
tron density scale height, and the core- mantle boundary in the Earth, which enters into 
the Earth effect calculation; the combined solutions were insensitive to the changes.0 
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Table 3. The best fit parameters of the combined MSW solutions. 





Nonadiabatic 
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3.3. Updated Results 

The updated MSW results including the most recent GALLEX, SAGE, and 
Kamiokande III data are presented in Fig. [7] (a). The joint fit includes the Homestake, 
combined gallium, Kamiokande II day-night data (6 data points) and the Kamiokande 
III data. We have incorporated the Earth effect, the theoretical and experimental un- 
certainties and their correlations. The Bahcall-Pinsonneault SSM with helium diffusion 
effect is assumed unless otherwise mentioned. The improved confidence level definition 
is robust under the assumption that the probability density is distributed throughout 
the parameter space and the distribution is described as a multiple gaussian distribu- 
tion, i.e., the allowed regions should have an elliptic shape, which is a good approxima- 
tion in our case. Fig. |7| (b) shows the combined allowed regions for 90, 95, and 99% C.L. 
Note that there is a third allowed region for Am 2 ~ 1.2 x 10~ 7 eV 2 and sin 2 29 ~ 0.8, 
but only at 99% C.L.; this region is disfavored by the Kamiokande II day-night data 
and also by the combined result of GALLEX and SAGE. 

The best fit parameters, x 2 values, the goodness-of-fit (G.O.F), and relative prob- 
abilities are listed for each solution in Table |3|. G.O.F is the probability of obtaining by 
chance a x 2 value equal to or larger than the obtained x 2 ■ The relative probability is 
the probability of finding the true solution in this region when the probability density 
is assumed to be a gaussian distribution for each solution.0 

Note that the small-angle (nonadiabatic) solution provides an excellent fit: the 
MSW solution is an acceptable theory by the combined observations. In fact x 2 — 2.4 
for 7 d.f. (9 data points - 2 parameters) is too good. This is due to a small contribution 
from the Kamiokande II 6 day-night data point (x 2 ~ 2.1), and without the day-night 
data, the x 2 A d.f. is 0.3 (58% C.L.), which is reasonable. 

Compared to the small-angle (nonadiabatic) solution, the two large-angle so- 
lutions are disfavored by the observations. This is due to the characteristic energy 
dependence of the solutions shown in Fig. For the small-angle (nonadiabatic) solu- 
tion, the 7 Be neutrinos and the lower energy part of the 8 B neutrinos are suppressed, 
consistent with the observed relative rates of Homestake and Kamiokande. On the 
other hand, the energy dependence for the large angle solution is mild, and the central 
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Fig. 7. (a) The MSW solutions allowed at 95% C.L. by the combined observations of the 
Homestake, the combined gallium SAGE and GALLEX, the Kamiokande II including six 
day-night data points, and the Kamiokande III data. Also shown are the regions allowed 
by each of the experiments and the region excluded by the Kamiokande II day-night data 
(dotted line), (b) The combined allowed regions at 90, 95, and 99% C.L. This is an updated 
analysis of Ref . pl| . 
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values of Homestake and Kamiokande cannot be simultaneously given by the MSW 
effect. This is similar to the reason why the astrophysical solutions are disfavored, but 
is less severe: in the large-angle MSW case, the converted or v T can still interact 
in the Kamiokande detector via the neutral current interaction (with about 1/6 the 
cross section), reducing the inferred v e survival probability to a value close to that of 
Homestake. There is also a slight enhancement of the survival probability at larger 
energies due to the Earth effect. 



o 




E v (MeV) 

Fig. 8. The energy dependence of the two MSW solutions allowed by the combined ob- 
servations. The small-angle (nonadiabatic) solution gives a better fit than the large-angle 
solution. The survival probability depends on the neutrino production distribution, and the 
plot assumes the distribution of the 8 B neutrinos, which is concentrated around the center 
of the Sun. The survival probability for the two 7 Be lines is essentially the same; for the pp 
neutrinos, the production region is more wide spread, and consequently the decrease of the 
survival probability around E ~ 0.5 MeV is less steep. 

The MSW oscillations also can be considered for sterile neutrinos, and the allowed 
parameter space is displayed in Fig. [|. There is no allowed large-angle region even at 
99% C.L., since the absence of the neutral current events in Kamiokande makes 
the v e survival probabilities effectively larger than for flavor oscillations. The large- 
angle regions are therefore strongly disfavored, much like the astrophysical solutions. 
The exclusion of the large-angle region for sterile neutrinos has been independently 
obtained from neutrino counting in big-bang nucleosynthesis.^ 

•t The MSW conversion formula is also modified for sterile neutrinos by the replacement n e — > n e —n n /2 
(see Eqn. ||), but the effect on the combined solutions are numerically negligible. 
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Fig. 9. The parameter space allowed for oscillations for sterile neutrinos. There is no allowed 
regions for large angles even at 99% C.L., which is consistent with the big-bang nucleosyn- 



thesis constraint. Taken from Ref. |31] 



The MSW regions for the Turck-Chieze-Lopes model is displayed in Fig. |10|. This 
model does not include the particle diffusion effect, and also the neutrino production 
distributions and the electron density distributions have not been published. The calcu- 
lation is done by using those of the Bahcall-Pinsonneauit model. The obtained allowed 
regions are slightly shifted outward of the MSW triangle compared to the Bahcall- 
Pinsonneauit solutions due to the smaller predicted fluxes. Also the allowed regions 
are noticeably larger, reflecting that the flux uncertainties are more conservative in the 
Turck-Chieze-Lopes SSM. Therefore the correct uncertainty estimation is important 
in constraining neutrino parameters, and underestimation of the flux uncertainties can 
lead to underestimation of the uncertainties of Am 2 and sin 2 29. 

3.4- MSW in Nonstandard Solar Models 

The SSM is assumed in the MSW analysis discussed in the previous section. 
Even if the MSW is proven, however, that does not confirm that the SSM is correct; 
the MSW effect can also take place in the nonstandard solar models. The SSM is still 
needs to be tested and calibrated experimentally. To address such a possibility, we have 
carried out the MSW analysis while the core temperature is allowed to change freely. 
The neutrino fluxes are changed according to the power law obtained by Bahcall and 
Ulrich (Eqn. []]), except that the T c dependence of the pp flux is obtained from the 
luminosity constraint (Eqn. Q). From a three parameter fit (Am 2 , sin 2 29, and Tc), we 
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Fig. 10. The allowed parameter space when the Turck-Chieze-Lopes SSM is assumed. Taken 
from Ref. 31]. 



constrained 

T c = 1.02 ±0.02 (la), (7) 

with Xmin = 2-1 for 6 degrees of freedom (d.f.), where Tc is in units of the Bahcall- 
Pinsonneault SSM (15.57 x 10 6 K). Interestingly the existing observations constrain 
the central temperature of the Sun at the 2% level even in the presence of the MSW 
effect. (Note that there is no consistent Tc to describe the observations without the 



MSW effect. See Section pTT] .) Also the value obtained for Tc is consistent with the 
Bahcall-Pinsonneault SSM value Tc = 1 ± 0.006. If the MSW is proven, the above will 
be the first observational constraint on the core temperature of the Sun at the 2% level. 
The combined MSW regions with free Tq are shown in Fig. [11]. 

Since the 8 B flux has the largest uncertainty (14% in the Bahcall-Pinsonneault 
SSM and 25% in the Turck-Chieze-Lopes model), we have also carried out the analysis 
with the 8 B flux as a free parameter. The constraint obtained from the observations is 
weak, but consistent with the SSM: 



0(B)/0(B) SSM = 1-43 + 0.65 - 0.42 (la). (8) 

The allowed regions for this case is shown in Fig. [T2|. 

We have also considered the MSW solutions for nonstandard solar models that are 
explicitly constructed, and the results are shown for four nonstandard solar models in 
Fig. [L3|. Two of them predict smaller initial fluxes compared to the SSM: the low opacity 
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Fig. 11. The allowed parameter space when the central temperature is allowed to change 
freely. The data constrain Tc = 1.02 ± 0.02 (lcr). Taken from Ref. |]31j]. 
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Fig. 12. The allowed parameter space when the 8 B flux is used as a free parameter. Taken 
from Ref. ||. 
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mo dell and the large Su model.ii The other two predict larger fluxes than the SSM: 
the high Y modelS and the maximum rate model.0 With nonstandard input parameters 
or nonstandard core temperatures, a wider MSW parameter space is allowed. 

Among the above nonstandard solar models, the maximum rate model and the 
large Su model are significantly different from the SSM, and are disfavored by the solar 
neutrino data even with the MSW effect. But nonstandard core temperature within 
Tq ~ 1 — 1.04 is perfectly consistent with the observation. A further constraint on Tq 
(or the neutrino flux) from the helioseismology observations is welcome. 

3. 5. Observational Constraints and Hints for Neutrino Mass and Mixing 

There are two observational hints of neutrino mass other than solar neutrinos: 
the oscillation interpretation of the atmospheric neutrino deficit ,0 and the cold plus 
hot dark matter scenario for the cosmic background radiation measurements and the 
large-scale structure observations.!^! The parameter ranges for those hypotheses are 
displayed in Fig. [14|, along with the MSW and vacuum oscillation solutions^ of the 
solar neutrino problem. The constraints from the oscillation experiments!! and the 
constraints for oscillations to sterile neutrinos from big-bang nucleosynthesis^ are also 
shown. 

Fig. |14| also shows the parameter space of various theoretical predictions involving 
the see-saw mechanisnfl the v e — mixing in the SOio grand unified theory (GUT) 
with an intermediate scale breaking,H@ the v e — v T mixing in the supersymmetric 
SO 10 GUT,B@ and the v e — mixing in the superstring-inspired model with non- 
renormalizable operators!! (there is no compelling prediction for the mixing angles for 
this model). The mixing angles shown for these theories (Fig. [T4|) are the corresponding 
quark mixing angles (Vi ept 

on ~ ^ckm)- In the SO\q GUT, the v T mass is expected to 
be in a range relevant to the cosmological hot dark matter. The detailed mass range 
and mixing angles are model dependent, and can be different from what is shown. 
The general idea of grand unification and the see-saw mechanism suggests the neutrino 
mass scales shown in Fig. ITJL and there is a solid motivation to look for neutrino mass 



in the MSW range relevant for solar neutrinos. 

The next question is whether this picture is reality, or just a fantasy of particle 
physicists, and, more importantly, whether we can prove or disprove such theoretical 
ideas. For solar neutrinos, experimental quests are directly aiming at the detection 
of neutrino oscillations, and the next generation experiments of the Sudbury Neutrino 
Observatory (SNO), Super-Kamiokande, BOREXINO, and ICARUS should either con- 
firm or rule out the neutrino oscillation hypothesis. For other neutrino mass ranges, 
two accelerator-based oscillation experiments at CERN (CHORUS and NOMAD) will 
start operating in the spring of 1994, and explore the cosmologically interesting v r 
mass range down to small mixing angles (Am 2 > 1 — 10 eV 2 and sin 2 28 < 2 x 10~ 4 
for <-> v T oscillations). The proposed long-baseline experiments of Brookhaven and 
Fermilab-SOUDAN will explore the parameter space suggested by the atmospheric 
neutrino deficit. 
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Fig. 13. The allowed parameter space for nonstandard solar models: (a) the low opacity 
model and (b) the large Su model predict smaller neutrino fluxes than the SSM, while (c) 
the high Y model and (d) the maximum rate model predict larger fluxes. Taken from Ref. pl|. 



25 




sin 2 29 



Fig. 14. Various observational constraints and hints for neutrino mass and mixing. Also 
shown are various theoretical predictions. See text for details. 



3.6. MSW Predictions for Future Experiments 

As discussed in the previous sections, we have compelling reasons to consider 
the MSW mechanism as a solution to the solar neutrino problem: a lack of viable 
astrophysical explanations provides a strong motivation for particle physics solutions, 
while the data are consistent with the MSW descriptions; the idea of GUTs and the 
see-saw mechanism suggests a neutrino mass range consistent with the MSW solutions. 

It is still premature, however, to convince ourselves that neutrino mass and mixing 
are the cause of the solar neutrino deficit. More data are needed, and the next genera- 
tion experiments such as SNofl Super-Kamiokande,il BOREXINO,0 and ICARUS0 
should provide a definitive answer for neutrino oscillation hypothesis. There are two 
theoretical questions which should be answered by the new experiments. 

• Calibrate the solar model. The neutrino flux should be determined component 
by component from experimental data. 

• Determine the solution and the parameter space. In particular distinguish the 
MSW solutions from astrophysical solutions and also from other particle physics 
solutions, and constrain Am 2 and sin 2 29. 

The procedure of doing these can be complicated, since the the determination of the 
neutrino parameters depends on the neutrino flux and, in turn, the determination 
of the neutrino flux requires a knowledge of the neutrino parameters. These can be 
accomplished by the next generation experiments with high counting rates and various 
reaction modes. The reaction modes of the experiments with projected counting rates 
are summarized in Table |j. 

Given the current data and the SSM, one can yield robust predictions for future 
experiments to test the MSW hypothesis. The MSW mechanism provides a solid, cal- 
culable theoretical frame work; the theoretical uncertainties are, by and large, under 
control. The predictions include the charged to neutral current ratio, the spectrum dis- 
tortions, and the day-night difference. The charged to neutral current ratio (CC/NC) 
is the gold-plated measurement in SNO of confirming neutrino oscillations; however, 
whether SNO can distinguish the MSW solution from astrophysical solutions at the 
3a level is solar model dependent; if the 8 B flux is smaller by 20% than that of the 
Bahcall-Pinsonneault SSM, the MSW solution may not be established given the pro- 
jected systematic uncertainty of ~ 20%. On the other hand, energy spectrum distortions 
and the day-night effect are completely free from the solar model uncertainties and the- 
oretically clean observables. By requiring consistency among those measurements from 
different experiments, we should be able to determine the solution and constrain the 
parameter space with a sufficient precision. 

3.6.1. CC/NC Measurement in SNO 

Measuring a depletion of the charged current (CC) rate with respect to the neutral 
current (NC) rate in SNO is considered to be the smoking gun evidence to establish 
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Table 4. The future solar neutrino experiments. CC, NC, and ES stand for charged current, 
neutral current, and elastic scattering mode. 



Experiment 


Commission date 


Reaction mode 


Event rate 


SNO 


Fall, 1995 


CC: v £ + d^ p H 


- p + e~ 


3,000 cvents/yr 






NC: i^, M ,T + d- 


■ ^e,n,T + P + n 


3,000 events/yr 






ES: u e ,n,T + e~ - 


~* ^e, M ,T + e~ 


200 events/yr 


Super-Kamiokandc 


April, 1996 


ES: Ue,n,r +e~ - 


-> ^e, M ,T + e~ 


8,000 events/yr 


BOREXINO 


Late 1990's 


ES ( 7 Be): u e ^, T 


+ e~ — > y eiMlT + e - 


18,000 events/yr (SSM) 


ICARUS 


~ 2000 


ES: Ve^.T + e~ - 


■+ ^e,,j,r + e~ 


2,900 events/yr /module 






CC: v e + 40 Ar -. 


- 40 K*+e- 


2,400 events/yr /module 


Iodine 


Late 1990's 


CC: u e + 127 I -» 


l27 Xe + e - 





neutrino oscillations. From the global solutions assuming the Bahcall-Pinsonneault 
SSM [Fig. (a)] at 95% C.L., the ratio (CC/NC) is expected to be 

(CC/NC) SNO f 0.2 - 0.65 (Nonadiabatic Solution) ( . 

(CC/NC) NoQsc ~ \ 0.2-0.3 (Large- Angle Solution), U 

where (CC/NC) No Osc is the ratio expected for no oscillations or for oscillations into 
sterile neutrinos. This prediction is, however, solar model dependent. Using the Turck- 
Chieze-Lopes model, which predicts a smaller 8 B flux with a larger uncertainty, the 
prediction [Fig. [15] (b)] becomes 

(CC/NC) SNO f 0.2-0.85 (Nonadiabatic Solution) 

(CC/NC) No Gsc ~ \ 0.2 - 0.4 (Large-Angle Solution). ( ' ' 

I emphasize that the prediction of whether SNO can distinguish the MSW solutions 
from astrophysical solutions by the CC/NC measurement is strongly dependent on the 
8 B flux and its uncertainty used in obtaining the current MSW solutions. With the 
current projection of systematic uncertainty at the 20% level, the distinction with a 
3a significance is achievable assuming the Bahcall-Pinsonneault SSM. It is no longer 
guaranteed if the Turck-Chieze-Lopes SSM is assumed. Given the large uncertainty in 
the 5*i7 determination, a 8 B flux smaller by 20% than that of the Bahcall-Pinsonneault 
SSM is a distinct possibility. Furthermore, there are additional uncertainties in other 
inputs such as opacity calculations. It is emphasized that the solar model uncertainties 
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can affect the predictions for future experiments, and especially the SNO experiment 
development can be affected by future improvements in the Sn determination. 




SNO (Charged Current) 

I I Combined Fit 95% C.L 



10 - 



10 



10"° 



> -6 

& 10 



10' 



10"° 



10 10 10 10 

sin 2 29 



10" J 



10 




SNO (Charged Current) 

I | Combined Fit 95% C.L. 

Turck-Chieze-Lopes SSM 



10 10 
sin 2 29 
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(b) Turck-Chieze-Lopes SSM 



Fig. 15. The contours of the signal-to-SSM ratio for the SNO charged current measurement 
with the combined allowed region assuming (a) the Bahcall-Pinsonneault SSM and (b) the 
Turck-Chieze-Lopes SSM. This is an updated result of Ref. 



3.6.2. Spectrum Distortions 

The 8 B spectrum measurements in SNO, Super-Kamiokande, and ICARUS provide 
a unique opportunity to look for energy dependent effects on low energy neutrinos 
with an unprecedented precision. Astrophysical effects on the spectrum such as from 
gravitational red-shifts and thermal fluctuations are completely negligible at the ob- 
servable level,0 while particle physics effects (or new neutrino physics) is often energy 
dependent, and spectrum distortions are a theoretically clean signal for new physics. In 
particular the MSW small-angle (nonadiabatic) solution predicts spectrum distortions, 
suppressing the lower energies of the 8 B flux more. The spectrum shapes predicted 
by the global MSW solutions are displayed in Fig. [16] for SNO (charged current re- 
action) and Super-Kamiokande. The error bars correspond to statistical uncertainties 
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from 6,000 and 16,000 events, respectively (equivalent to 2 years of operation). To com- 
pare differences in the shape, the large-angle spectrum is normalized to the small-angle 
(nonadiabatic) spectrum above the threshold (5 MeV). The calculation includes the 
detector energy resolutions and also the differential charged current cross section@@ 
for SNO. [j] The small-angle (nonadiabatic) solution shows a depletion of the spectrum 
at lower energies (< 7 MeV), and it is essential to achieve a low threshold to observe the 
distortion. The large-angle solution shows little distortion, and the shape is identical 
to astrophysical solutions (or no oscillation) in the scale of the figure. 



3.6.3. Day-Night Effect 

The day-night effectffl is another purely particle physics mechanism that is free from 
astrophysical uncertainties since the day-night difference (or time dependence) of the 
signal is independent of the absolute 8 B flux. The regeneration of u e 's in the Earth can 
be measured not only by the day-night difference, but also by time-dependence during 
the night and energy spectrum distortions. The time dependence and the spectrum 
distortions are often strongly correlated. The parameter space for which the day-night 
effect is most drastic has already been excluded by the Kamiokande II data, but the 
next generation experiments with higher statistics are sensitive to a wider parameter 
space. The day-night effect, if observed, is a dream observable for particle physicists: 
the determination of Am 2 and sin 2 29 can be done precisely by measuring the time 
dependence and spectrum distortions, completely independent of solar model predic- 
tions. The calibrations of the solar neutrino fluxes could then be easily done once Am 2 
and sin 2 29 are determined from the day-night effect. 

The large-angle solution of the combined analysis predicts a night-time enhance- 
ment of the signal which should be observable in SNO and Super-Kamiokande. The 
predictions of the day-night effect in the small-angle (nonadiabatic) and the large-angle 



solutions are shown in Fig. |17| and [Tj|. In fact the enhancement is such a large effect 
that it might be observable in the combined data of Kamiokande II and III (Fig. p"9|). 

An interesting possibility of detecting the day-night effect for the small-angle 
(nonadiabatic) solution was carefully investigated by Baltz and Weneser recently." It 
was pointed out that, for the best fit parameters in the small-angle (nonadiabatic) 
region, the 8 B neutrinos can resonate in the core of the Earth, which corresponds to 
the last bin in Figs. |17] (a) and [l^ (a). Separating those signals, it was concluded 
that Super-Kamiokande should be able to see the Earth effect at about 2a in a year. 
It is still an open question whether additional information from the energy spectrum 
distortions enhances the signals, and whether the effect is therefore observable with 
definitive statistics above the background. 



3 It is important in the spectrum shape analysis to include the detector energy resolutions and the 
charged current electron spectrum. The approximation "electron energy = neutrino energy - 1.44 
MeV" significantly overestimates the sensitivity to distortions. 
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SNO CC Electron Spectrum 
6000 Events (2 year) 



■ Nonadiabatic 

O Large-Angle / No Oscillations 




E e (MeV) 



(a) SNO (CC) 




(b) Super-Kamiokande 



Fig. 16. The spectrum shape expected for (a) SNO and (b) Super-Kamiokande. The error 
bars show the statistical uncertainties from totals of 6,000 and 16,000 events, respectively 
(equivalent to 2 years of operation). The large-angle spectrum is normalized to the nonadia- 
batic spectrum above 5 MeV. Taken from Ref. [31]. 
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Fig. 17. The day- night effect in SNO for (a) the small- angle (nonadiabatic) solution and (b) 
the large-angle solution. 5 sun is the angle between the nadir at the detector and the direction 
in the Sun. The error bars are the statistical uncertainties assuming a total of 3,000 events 



(equivalent to one year of operation). Taken from Ref. [31]. 
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Fig. 18. Same as Fig. 17, but for Super-Kamiokande. The error bars are the statistical 
uncertainties assuming a total of 8,000 events (equivalent to one year of operation). Taken 



from Ref. 1 31 ] . 
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Kamiokande Day-Night Data 
Large-Angle Solution 




Fig. 19. Same as Fig. 17, but for Kamiokande. The error bars are the statistical uncertainties 



assuming a total of 200 events. Taken from Ref. [31| 
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Fig. 20. The contours of the signal-to-SSM ratio for (a) Super-Kamiokande and (b) BOREX- 
INO, along with the global MSW solutions (from present data) obtained in the Bahcall- 
Pinsonneault SSM. This is an updated result of Ref. pi] . 
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3.6.4. Summary of MSW Signals 



The parameter regions allowed by the combined observation are shown with the con- 
tours of the signal-to-SSM ratio in Fig. |20| for Super-Kamiokande and BOREXINO, 
assuming the Bahcall-Pinsonneault SSM. Given the high event rates in SNO, Super- 
Kamiokande, and BOREXINO, the accurate determination of the neutrino parameters 
should be possible. 

The various characteristic signals from the global MSW solutions are summarized 
in Table [5| Also listed is the vacuum oscillation solution (Am 2 ~ 8 x 10~ n eV 2 and 
sin 2 29 ~ 0.8), which is allowed by the data, but requires a fine-tuning between the 
Sun-Earth distance and the neutrino parameters. It is stressed that the significance of 
establishing neutrino oscillations by the CC/NC measurement in SNO depends on the 
solar model. On the other hand, the effects of spectrum shape distortions, the day-night 
effect, and seasonal variations are free from solar model uncertainties and theoretically 
clean. An equivalent list for oscillations to sterile neutrinos are shown in Table IS. 



Table 5. Characteristic signals expected in the global MSW solutions in flavor oscillations. 
The second large-angle solution is strongly disfavored by the data, and the solution only ap- 
pears at 99% C.L. The vacuum oscillation solution is also listed. Whether SNO can definitively 
establish neutrino oscillations by the CC/NC measurement is solar model dependent. 
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Table 6. Characteristic signals expected in the global MSW solutions for sterile neutrinos. 
The large- angle solutions (both (I) and (II)) and the vacuum solutions are strongly disfavored 
by the data and excluded for sterile neutrinos. The large-angle (I) solution is independently 
excluded by the big-bang nucleosynthesis constraint 
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4. Summary and Wish List 

4-1. Summary 

Various theoretical models were confronted with the solar neutrino observations as 
possible solutions of the solar neutrino deficit. The astrophysical solutions are strongly 
disfavored by the data: (a) The SSM is excluded by the each of the experiments, (b) 
The model independent analysis concludes that the combined results from Homes- 
take and Kamiokande exclude any astrophysical solutions, (c) Even if the Homestake 
data were discarded, we have no realistic solar model that simultaneously explains the 
Kamiokande and the gallium results within the experimental uncertainties. 

The MSW mechanism, on the other hand, provides a complete description of the 
data, without discarding any of the observations. The parameter space allowed by the 
data is constrained to 

Nonadiabatic: Am 2 ~ 10~ 5 eV 2 and sin 2 29 ~ 0.007 (11) 
Large-angle: Am 2 ~ 10" 5 eV 2 and sin 2 29 ~ 0.8. (12) 

There is a third allowed region for Am 2 ~ 10~ 7 eV 2 and sin 2 29 ~ 0.8, but only at 99% 
C.L. For oscillations to sterile neutrinos, the solution is allowed only in the nonadiabatic 
region. 

The obtained parameters are sensitive to the solar model uncertainties. When 
nonstandard core temperatures are allowed, the data constrain Tq = 1.02 ± 0.02 (la), 
consistent with the SSM. A wider MSW parameter space is possible in nonstandard 
solar models. 

The predictions for the CC/NC ratio, spectrum distortions, and day- night effect 
in SNO and Super-Kamiokande were discussed in detail. It is stressed that the MSW 
prediction for the CC/NC ratio is solar model dependent, and a smaller 8 B flux can 
make the difference between the MSW solutions and no-oscillations ambiguous. 
4.2. Wish List 

The solar neutrino experiments have provided us an opportunity for a better 
understanding of the Sun as well as a possible discovery of new neutrino physics. There 
will be more data soon, and much work is still needed involving experimentalists, 
astrophysicists, nuclear physicists, as well as particle physicists. The INT workshop 
allowed us to convene and discuss this inter-disciplinary topic, and, taking advantage 
of the opportunity, I will address my hopes from a particle physicist point of view: 

From experimentalist friends: 

• Uncertainty correlation matrix be published. In constraining the parameters in the 
MSW analysis, the energy spectrum and the time-divided data are very informa- 
tive. From published papers, however, it is often not clear which of the systematic 
uncertainties are correlated among data points. In high-energy experiments such 
as in LEP, it is becoming more common to publish the correlation matrix, and I 
hope this trend will be followed in solar neutrino experiments. 
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Prom astrophysicist friends: 



The Syi uncertainty be reevaluated. It is my belief that the current uncertainty ob- 



tained in Ref. [15j and used in the Bahcall-Pinsonneault model is underestimated. 
A larger uncertainty does not solve the solar neutrino problem; however, it does 
change the obtained MSW parameter space and affects the MSW predictions for 
the future experiments, especially the CC/NC ratio in SNO. 

• A search for nonstandard solar models consistent with the Kamiokande and gal- 
lium data be pursued. There is no way one can reconcile the Homestake and 
Kamiokande data by astrophysical explanations. Even discarding the Homes- 
take results entirely, we have no realistic solar model consistent both with the 
Kamiokande and gallium data, but more studies are needed. Especially, every 
nonstandard solar model should be confronted with the helioseismology data. If 
we cannot find a solar model consistent with the Kamiokande, gallium, and the 
helioseismology data after discarding the Homestake results, our motivation for 
new neutrino physics will be complete. 

• Can the neutrino flux and the MSW parameters be constrained simultaneously? 
Once neutrino oscillations are established in the next generation experiments, 
the most significant task is to constrain the neutrino flux and the MSW parame- 
ters simultaneously from observations. This can be complicated since the MSW 
mechanism can distort the neutrino spectrum. Either simple relations among the 
fluxes (like the Tc power laws) or Monte-Carlo models like the ones by Bahcall 
and Ulrich but with a wider range of solar models are needed. It would be ideal 
if the helioseismology data were incorporated in such an analysis. 
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